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ABSTRACT 

Adaptive optics (AO) systems have significantly improved astronomical imaging capabilities over 
the last decade, and are revolutionizing the kinds of science possible with 4-5m class ground-based 
telescopes. A thorough understanding of AO system performance at the telescope can enable new 
frontiers of science as observations push AO systems to their performance limits. We look at recent 
advances with wave front reconstruction (WFR) on the Advanced Electro-Optical System (AEOS) 
3.6 m telescope to show how progress made in improving WFR can be measured directly in improved 
science images. We describe how a "waffle mode" wave front error (which is not sensed by a Fried 
geometry Shack-Hartmann wave front sensor) affects the AO point-spread function (PSF). We model 
details of AEOS AO to simulate a PSF which matches the actual AO PSF in the /-band, and show 
that while the older observed AEOS PSF contained several times more waffle error than expected, 
improved WFR techniques noticeably improve AEOS AO performance. We estimate the impact of 
these improved WFRs on if-band imaging at AEOS, chosen based on the optimization of the Lyot 
Project near-infrared coronagraph at this bandpass. 

Subject headings: instrumentation: adaptive optics - instrumentation: miscellaneous - techniques: 
miscellaneous - methods: numerical 



1. INTRODUCTION 

Various wave front sensors (WFS's) are in principle 
insensitive to certain kinds of wave front errors. The 
comm only- used square 'Fried geometry' sensing (jFriedl 
Il977|) in most Shack-Hartmann-based adaptive optics 
(AO) systems is insensitive to a checkerboard- like pat- 
tern of phase error called waffle. This zero mean slope 
phase error is low over one set of WFS subapertures ar- 
ranged as the black squares of a checkerboard, and high 
over the white squares. 

Since waffle modes are in the null space of a Fried ge- 
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ometry WFS, they can be present in the reconstructed 
wave front unless the wave front reconstructor (WFR) 
explicitly removes or attenuates these modes. Because 
waffle is so common to Shack-Hartmann systems - 
early images from th e 241-actuator Palomar AO sys- 
tem llTrov et al.ll2000ft the 349 -actuator Keck AO sys- 
tem ifv^zrnowich et alJ teOOObj) . and the 941-actuator 
Advanced Electro Optica l System (AEOS) AO system 
ijRoberts fc Nevmanl2002|) all showed the presence of sig- 
nificant waffle — effort has been made to understand how 
to control this and other "blind" reconstructor modes 
lIGavell 120031 IPovneerl 120031 IMakidon et alJ 12003). Ad- 
vances in wave front reconstruction approaches have im- 
proved AO system performance considerably on all three 
of the aforementioned telescopes (M. Troy 2003, private 
communication; D. Gavel 2003, private communication; 
G. Smith 2003, private communication). 

However, the prevalence of waffle in most Shack- 
Hartmann-based AO systems led us to investigate the 
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physical nature of waffle, without a sophisticated control- 
theoretic formalism. Here we present the understanding 
we gained about waffle mode error using a combination 
of Fourier analysis and numerical simulation to examine 
the effects of waffle on the observed point spread func- 
tion (PSF). We concentrate our efforts on the fundamen- 
tal waffle mode, where the checkerboard extends over the 
full aperture of the telescope (a c omplete enum eration of 
all waffle modes can be found i nlOavel ll200ll. 

We follow the methods of iSivaramakrishnan et alJ 
1)20011) (hereafter referred to as SKMBK), which estab- 
lished simulations calibrated to AO imaging data from 
Palomar as a viable means to examine AO system char- 
acteristics. We have tuned our simulations to the AEOS 
Adaptive Optics PSF in an attempt to reproduce waffle 
behavior seen in early data acquired with this system. 

The 3.6 m AEOS telescope, part of the U.S. Air Force 
Research Laboratory's Maui Space Surveillance System 
(MSSS), is arguably the best telescope on which to 
develop the nascent field of Extreme Adaptive Optics 
(ExAO). This is a result both of the telescope's location 
at a prime astronomical site possessing good seeing and, 
more importantly, the AEOS adaptive optics system's 
941-actuator deformable mirr or, with 35 actuators span - 
ning the diameter of the DM ({Roberts fc Nevmarl 2002'). 
This provides the highest actuator density available to 
any civilian astronomical observing program, with a pro- 
jected actuator spacing of order 0.11 m/actuator at the 
primary mirror. 

In Section |3 we examine the nature of waffle mode 
error using Fourier ana l ysis a n d the PSF expansi on of 
ISivaramakrishnan et all l)2002|k iPerrin et all ( |2003|) . We 
apply this understanding to simulations of waffle mode in 
AEOS /-band images in Section 13 We then extrapolate 
our /-band imaging simulations to iZ-band in Section 0] 
to predict AO system performance in the near-infrared 
at AEOS. Our choice in simulating //-band is predicated 
on the Lyot Project coronagraphic imager's optimization 
for operations at JJ-band and the scientifi c potential of 
this instrument l)Onpenheimer et al.ll2004L 4 

2. THE NATURE OF WAFFLE MODE ERROR 

Waffle mode error is inherent in any wave front sensing 
that uses a square array of wave front slope sensors. The 
atmospheric phase disturbance typically possesses some 
power in this mode, which will "leak" through the AO 
system to produce a square grid of ghost images centered 
around the true PSF peak. We examine this process in 
detail, first for a simplified case consisting of a single 
WFS cell and then extending to the full sensor array. 

The Shack-Hartmann WFS geometry is the most com- 
mon choice for AO system with more than a few tens of 
actuators. Its choice is driven by the essentially square 
nature of WFS detector pixel layout, as well as the rel- 
ative simplicity of th e concept. An other type of WFS, 
the curvature sensor fRad-difij 11988), is not sensitive to 
waffle mode error, although it has its own particular null 
space. Curvature sensors are typically not used for AO 
systems with more sensin g channels because of thei r noise 
propagation prop erties llGlind emann et alJ 120 00). AO 
syste ms at Keck jW izinow ich et alJ l2000albfl . Palomar 
Hale i|Trov et aLlfeOPOj) . and AEOS ^Roberts k Nevmanl 

4 See also: http://www.lyot.org 



l2002h all use Shack-Hartmann WFSs, imaging light onto 
a CCD detector using square lenslet arrays located at 
an image of the entrance pupil (the Shane 3 m telescope 
at Lick Observatory is an exception, having a hexagonal 
grid Shack-Hartmann WFS geometry). The WFS de- 
tector is read out, and the on-chip x and y coordinates 
of the spots' centroids are delivered to the wave front 
reconstruction computer. The difference of one of the 
centroid's coordinates from its long-term temporal mean 
value is directly proportional to the slope s = [s x ,s y ) 
of the wave front over the corresponding subaperture. 
The constant of proportionality is just the plate scale 
of the WFS lenslet. We can therefore use the quantity 
s = (s x ,s y ) = (n(x c - (x c )),K,(y c - (j/c))) as the mea- 
sured quantity in any given subaperture. Here (x c ,y c ) 
is the location of an instantaneous centroid of a spot in 
the Shack-Hartmann subaperture (in pixels), () denotes 
temporal averaging, and k is the plate scale of the detec- 
tor (in arcseconds or microradians per pixel). This gives 
us the vector slope of the wave front (in angular units) 
on a grid of locations over the full aperture. 

2.1. The unit cell 

The nature of the smallest spatial scale waffle mode 
can be deduced from inspecting a 'unit cell' of the WFS 
array. We assume a WFS subaperture which has four 
DM actuators, one at each corner of a square (see Fig- 
ure ^) . The relationship between the wave front slope 
across this subaperture and the heights of wave front at 
the four corners of the square contains the key to the waf- 
fle mode error, since the DM actuators at these corners 
will need to be moved by exactly half these heights to 
flatten the wave front across the subaperture. Of course 
the actuators could be placed anywhere relative to the 
subaperture grid, though real AO systems typically try 
to maintain the DM actuator-to-pupil alignment we de- 
scribe if the actuator spacing and geometry matches the 
subaperture geometry. We stress here that waffle mode 
error does not depend fundamentally on the actuator lo- 
cation and geometry: it only depends on the sensor ge- 
ometry. 

We label the corners of the subaperture with mnemon- 
ics generated from the cardinal directions, assuming 
North points toward the top of the figure. The corner 
heights are therefore h = {hNE,hNW,hsE,hsw)- We 
assume that the wave front is flat across the subaper- 
ture, so we can immediately deduce the way our slope 
measurement s constrains the wave front heights at the 
corners, 

(h NE + h S E)/^d- (h NW + hsw)/2d=s x 
(h NW + h NE )/2d- (h S E + hsw)/2d=s y , 

where d is the subaperture spacing as well as the inter- 
actuator spacing. These are the two equations of con- 
dition of the problem. This relation between the four 
unknowns h and the two measured quantities s can be 
re-expressed as 

hsE — hNW = d(s x — s y ) 
hNE - h S w = d(s x + s y ). 

We use our only two observable quantities, s x and s y , 
to relate the heights across the diagonals of the subaper- 
ture. This leaves two undetermined quantities in these 
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equations relating our unknown heights to the two data 
points: the equations of condition are rank-deficient by 2. 

It is not possible to determine all four heights in h 
with just two observables. If we face the engineering 
problem of having to command four actuators based on 
the data contained in s, we will have to add two arbitrary 
equations to the defining equations of the problem to 
solve them. That is, we have to add two more equations 
to the system of equations in order to obtain a solution 
for h. One such equation is easy to come up with: the 
physics of the problem indicates that the mean value of 
the four heights can be any value we choose, as it is 
the piston of the wave front. We can therefore add the 
constraint equation 

(h NE + h NW + h S E + hsw)/^ = P (1) 

to the equations of condition. Here P, the mean height 
of the wave front, is arbitrary. 

Our original defining equations are difference equa- 
tions, so it is scarcely surprising that the offset of the 
heights is unknown. Wave front slope data alone cannot 
determine the piston of the wave front. Let us assume 
a zero piston in our case (we will set the mean voltages 
sent to the DM actuators using common sense in any real 
system) : 

hNE + h N w + hsE + hsw — 0. (2) 

We still need one more equation without data — a 
constraint equation — to solve our problem. This is the 
equation that determines how much waffle we introduce 
into our very simple AO system. So far we have set the 
slopes across the unit cell diagonals using data, and the 
mean piston using common-sense. We must now guess 
how the means of the two sets of diagonally-related actu- 
ators differ from each other, in the absence of new infor- 
mation (data). If the wave front is relatively flat across 
the sub aperture (e .g., d <C rg, where rg is the Fried 
length l)Friedl H9661). we can assert that the means are 
equal. This statement, expressed as the equation 

(h N E + h S w) ~ {hNW + hss) = 0, (3) 

adds the equation required to solve the system of equa- 
tions relating our two WFS slopes to our four actuator 
heights completely. It determines how much waffle mode 
present in the atmosphere we let slip through our AO 
system, on average. Equation Q is an equation of con- 
straint, since it does not involve any data. We note that 
ellipticity of the individual WFS subapertures's spots can 
yield data on the waffle mode error, although no AO sys- 
tem we know of uses this information in practice. 

A singular value decomposition (SVD) method is of- 
ten used to solve the underdetermined set of equations 
between WFS data and actuator heights. The physics 
of what really happens when we perform a wave front 
reconstruction with SVD is not apparent: this method 
does not distinguish between constraints and equations 
of condition. 

2.2. The full lattice 

Here we extend the description in Section 12.11 to the 
simplest case when we have more data than unknowns. 
We will use a square aperture. While this is obviously 
unrealistic, our analysis of this case can be applied to 



arbitrary aperture geometries, and leads to an under- 
standing of the cause of unsensed modes in wave front 
reconstructors that use slope measurements as input to 
their algorithms. 

Any Fried geometry WFS system on a simply- 
connected aperture will possess two families of actuator 
heights. Extending the unit cell to a case when there are 
more data points than unknown actuator heights demon- 
strates this. The N s WFS subapertures produce 2N S 
slope measurements, since each subaperture yields an x 
and a y slope. With a square aperture, the simplest such 
case is that of a grid of four actuators on a side, enclosing 
a square of subapertures three to a side (see Figure Q). 
We must deduce 16 actuator heights using 18 pieces of 
slope data. On the face of it the problem appears to 
be overdetermined, but is in fact still underdetermined. 
Writing out the equations that relate slopes to actua- 
tor heights explicitly, the set of equations governing one 
family of heights is 
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and the equations for the second family are 
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where 

S{ di^Siy Six) 

<Ji = d(siy + s ix ). (6) 

These two sets of nine equations each are disjoint - 
no individual unknown actuator height appears in both 
sets. We see by inspection that we can choose any one of 
the first set of heights arbitrarily. Setting hi to a fixed 
number sets every actuator height in the first set. Simi- 
larly, we can choose h 2 to be any value we please. Each 
of the two families of actuator heights therefore has an 
arbitrary net piston of its own (i.e., the mean of all 
the actuator heights in the family) . This is the common- 
sense expression of the statement that the matrix de- 
scribing the system is rank-deficient by 2. These pistons 
can individually be set to be equal, thereby selecting a 
'waffle-free' solution. This solution still has arbitrary to- 
tal piston. Again, we can select a zero piston solution 
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from all the waffle-free solutions by requiring the pistons 
of each family be zero: 



16 



0. 



(7) 



0=1 



The atmospheric phase disturbance at any given in- 
stant does not typically possess zero waffle when viewed 
through the WFS subapertures. That is, if one regards 
the square WFS lenslets (projected back to the entrance 
pupil) as a checkerboard, the mean phase over the black 
squares will not equa l the mean p hase over the white 
squares. However, as iGavell ()200rjf) shows, for subaper- 
tures that are small compared to ro the net waffle will 
be small compared to unity (waffle mode error has units 
of radians, since it is a mean of a phase). By enforcing a 
zero- waffle solution on our actuator heights we introduce 
some error in our fit of the wave front. This error pro- 
duces the familiar waffle pattern appearing at an angular 
radius of v / 2(A/2d) in the image (see Section 0J. 

2.3. Waffle Mode Error and the PSF 

Here we consider waffle mode error as a phase aber- 
ration 4> w (x,y) over the pupil, and trace the effect 
of the waffle mode on the AO-corrected monochro- 
matic wave front analy t ically. We use the notation of 
iSivaramakrishnan et alJ l)2002jh which we restate here. 

The telescope entrance aperture and all phase effects 
in a monochromatic wave front impinging on an opti- 
cal system can be described by a real aperture illumina- 
tion function A(x,y) multiplied by a unit modulus func- 
tion A<f,(x, y) — e l ^ x - y \ The aperture plane coordinates 
(x, y) are in units of the wavelength of the monochro- 
matic light under consideration (to simplify the Fourier 
transforms), while the corresponding image plane co- 
ordinates (£,77) are in radians. Here, phase variations 
induced by the atmosphere or imperfect optics are de- 
scribed by a real wave front phase function (j) which pos- 
sesses a zero mean value over the aperture plane, 



J A(j)dxdy 
J Adxdy 



0. 



(8) 



We chose the location of the image plane origin at the 
centroid of the image PSF, which corresponds t o a zero 
mean tilt of the wave front over t he aperture ijTeagud 
119821: ISivaramakrishnan et al.lll995j) . 

We assume that the transverse electric field in the im- 
age plane is given by the Fourier transform (FT) of the 
field in the aperture plane l)Goodmanlll995|) . We indicate 
transform pairs by a change of case, where the FT of A 
is a and the FT of <fr is $. 

The "amplitude spread function" (ASF) of an optical 
system with phase aberrations is the FT of A^ = AA^, 
and is given by — a * a^, where * denotes convo- 
lution. This complex-valued quantity is proportional to 
the amplitude of the incoming wave in the image plane. 
The complex number's phase at any point is interpreted 
as a relative phase difference in the wave front at that 
point in the image plane. The corresponding PSF of this 
optical system is p^ = a^aj, where the superscripted * 
indicates complex conjugation (see e. o. . iBracewelll (2000) 
for the fundamental Fourier theory and Fourier optics 
conventions). 



For this analysis, we assume that the total phase 
<j>(x,y) is given by 4>ao (x,y)+(f> w (x, y), where (j>Ao(x,y) 
is the residual atmospheric phase after correction by the 
AO system. Thus the unit modulus function A^x^y) 
can be written as 

A = AaoA^ = e ^o(z,y)+<?Uz,y)]. (9) 

In order to understand the effects of waffle mode error on 
the PSF, we now look at a waffle aberration without any 
atmospheric aberration contribution. Setting 4>ao = 0, 
the aperture illumination function reduces to 

A$ = AA^^ = Ae i ^ x ' v \ (10) 

In Section 12.21 we presented waffle mode error as a 
checkerboard pattern in phase with a period of 2d, where 
d is the WFS subaperture size in the Fried geometry. In 
the one-dimensional case, we write this as a top-hat func- 
tion hU(x/d) convolved with series of Dirac <5-functions 
spaced at a period 2c?. Here h is a function of the atmo- 
spheric phase 4> a tm , defining the "strength" of the waffle 
function. The periodic function describing our checker- 
board thus takes the form 

n 

where we use the shah function, III(x/2d) to denote 
a series of Dira c S- functions at the required spacing 
ijBrace wcll 200(J. The waffle phase error (f> w can then 
be written 



nf- 



III 



(2d) 



* hIL 



(12) 



where D is the telescope aperture diameter in units of 
wavelength. The wave amplitude in the image plane due 
to waffle mode is the FT of (f> w : 

<5> w = hD smc(D8) * [2d III (2d0) • d sinc(d0)] . (13) 

The above expression is the field due to the telescope 
aperture convolved with a sine function sampled at 
8 = n/2d, where n is an integer. However, the func- 
tion sinc(c£0) is zero for each even value of n (e.g., 
n = 2, 4, 6, . . .). Thus this fundamental waffle mode error 
results in the replication of the perfect PSF at a lattice of 
points at the angular period equal to the spatial Nyquist 
frequency of the WFS subaperture density, the relative 
intensity of each point being modulated by sva.c(d6) (with 
half of those points falling on the zeros of the sinc(d#)). 

2.3.1. The PSF Expansion 

iPerrin e~ al. (2003) developed a general form for an 
expansion of an aberrated PSF, p, in terms of the Fourier 
transform of the phase aberration <f> over the aperture. 
They express the PSF was given as a sum of individual 
terms in an infinite convergent series, with the n th term 
of that series given by 



P 



(-1) 



n — k 



k=0 



k\(n-k)\ 



(a ** $)(a** n - fc $*), (14) 



where * n is used to denote an n-fold convolution op- 
erator. We refer the reader to the discussion of the 
se cond-order expansion of t he pa rtially co rrected PSF 
in ISivaramakrishnan et all l)2002|) and in IPerrin et alJ 
(2003), noting here that while the first order term 



Pi 



a*(a * $) - a(a* * $*) 



(15) 
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is due entirely to the antisymmetric component of A^, 
it does not contribute significantly to the waffle pat- 
tern. This is because the function a(x,y) is very small at 
points beyond a few A/ D from the center of the on-axis 
PSF, and the f i rst or d er term is multiplied by a (see 
iBloemhof et all (|200lD : iSivaramakrishnan et all ( 2002) 
for further discussion of pinned speckles) . This is borne 
out by simulations with pure waffle-mode phase aberra- 
tions showing the symmetric term s of the PSF expan- 
sion (e.g., the P2 and p^ terms of iPerrin et alJ ([2003)) 
dominate, even when the phase error itself is large and 
antisymmetric. We discuss this further in Section 13.1.21 

3. SIMULATIONS 

Our goal in this analysis is to understand the current 
performance of the AEOS AO system in the /-band, and 
to enable predictions of AO system performance in the 
-ff-band. Such predictions are useful for determining 
the contrast ratios accessible to the Lyot Project AO- 
optimized near-infrared (NIR) coronagraph on AEOS. 
Our simulations incorporate the AEOS pupil geometry, 
namely a 3.63 m clear aperture for the primary mirror, 
with a 0.40 m obstruction from the secondary and ter- 
tiary mirrors. We do not model secondary support struc- 
tures. 

A detailed simulation of an AO system including all 
aspects of wave front sensing, reconstruction, and cor- 
rection requires substantial computational effort. For 
some purposes it is sufficient to model the action of an 
AO system as a transfer function or spatial frequency 
filter acting on the incident wave front. In SKMBK, we 
used a high pass filter in spatial frequency space for this 
purpose. We tuned our initial mo del to fit data from 
the Pa lomar AO system reported in lQppenheimer et alJ 
( 2000). We briefly describe the SKMBK model here, re- 
peating some of the arguments presented in that work, 
and show how we modified the model to fit the AO imag- 
ing data from AEOS. 

The cutoff spatial frequency of the AO filter cannot be 
higher than the spatial Nyquist frequency of the actuator 
spacing, which is 

k AO = N act /2D = l/2d, (16) 

where N act = D/d, with D the primary mirror diameter 
and d the actuator spacing projected onto the primary 
mirror. At a wavelength A, this spacing in the pupil plane 
corresponds to a limiting angle of correction in the image 
plane, 

0ao = N act X/2D, (17) 

which defines the 'AO control r adius' ijPerrin et alJ 
(2003); P ovneer fc Macintosh! l|2004D discuss the AO con- 
trol radius in more detail). 

As we note in SKMBK, the shape of the AO filter near 
the cutoff frequency depends on details of the AO system. 
In general, DM actuator influence functions which extend 
to neighboring actuators positions reduce the sharpness 
of the cutoff. Noisy wave front sensing and intrinsic pho- 
ton noise reduce the efficacy of high spatial frequency 
wave front correction. The flow of the atmosphere past 
the telescope pupil and a non-isotropic refractive index 
spatio-temporal distribution also change the shape of the 
AO filter, as does imperfect DM calibration. As a result, 
the observed AO control radius is often smaller (in angle) 
than the theoretically- predicted AO control radius. 



Following the methods in SKMBK, we use a high-pass 
filter with a power law distribution in spatial frequency 
space to mimic the action of the AEOS AO system, 

p r k \ _ f {k/k AO ) n , ioik<k AO / lg N 
' 1 1, otherwise. 

In previous work ijMakidon et al.ir2003|) . we assumed a 
'parabolic' filter (e.g., n = 2) to model this AO sys- 
tem. However, further refinement of our simulations 
have shown us that a n = 0.9 power law distribution 
best fits the AEOS system (see Section 13.211 . We as- 
sume single-layer atmospheric turbulence, and generate 
independent realizations of Kolmogorov-spectrum atmo- 
spheric phase screens in 1024 x 1024 arrays to simu- 
late this turbulence. We use code implementing a rapid 
Markov metho d for gene rating Ko lmogorov-spectrum 
phase screens l |Lane et al.lll992t iGlindemann fc Daintvl 
119931 iPorro et al.ll2000ri ~with a spatial scale of 141.047 
pixels per meter in the aperture plane. 

We Fourier transform the input phase arrays, and mul- 
tiply them by the power law filter in Equation (|18l) to 
mimic the action of AO. Then we reverse-transform the 
spatially filtered arrays to obtain the AO-corrected wave 
front without waffle error. We use only the central half 
of the smoothed phase screens, in order to avoid any 
edge effects introduced by the finite extent of the phase 
screens during the numerical manipulations. 

3.1. Simulating Waffle Mode Phase Error 

For traditional Shack-Hartmann AO systems and 
WFR, waffle mode phase error present in the atmo- 
sphere remains unsensed and uncorrected, and thus leaks 
through the system. This is also true of the spatial filter 
we use to represent our AO correction. However, because 
of our code structure, we found it easier to add the waf- 
fle mode phase error back into the AO-corrected residual 
phase after correction rather than before. We found no 
difference between these two approaches. 

We calculated the waffle mode error, <fi w , present in 
our incoming atmospheric phase aberrations <j) a tm by first 
representing the AEOS WFS subapertures as a checker- 
board pattern over our aperture, with white and black 
squares corresponding to the two families of DM actua- 
tors (see Figure 0J- We then calculated the mean phase 
of the incoming wave front 4> a tm over the white squares 
(4>i ), as well as the mean over the black squares (4>2)- The 
zero-mean waffle phase error present in the atmospheric 
phase aberration is then constructed by assigning a phase 
4> w defined by (<f>i + <j>i)/2 over the white squares, and 
(<p\ — 4>i)/2 over the black squares. We then added <f> w 
back into our AO-corrected phase arrays (4>ao ) , produc- 
ing an AO-corrected wavefront with waffle mode phase 
error. 

Since 4>ao is nonzero in this case, the aperture illumi- 
nation function takes the form 

A^ = AA^ =J 4 e i [^( I .v)+*»(^!')]. (19) 

The PSF corresponding to this illumination function is 
just the square of the absolute value of the Fourier Trans- 
form of a^. 

We show this process graphically in Figure 01 over a 
portion of the AEOS pupil. Wc calculate the amount 
of zero-mean waffle phase error present in the incoming 
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phase screen (left) over the checkerboard array of sub- 
apertures. Without waffle mode phase error, the AO- 
corrected wave front would just be dominated by residu- 
als with spatial frequencies beyond those correctable by 
our spatial filter. However, the addition of waffle mode 
phase error produces a noticeable checkerboard pattern 
in the final, AO-corrected wave front (right). 

3.1.1. Waffle strength a function of tq 

As part of our analysis, we examined how <f> w varies 
with the Fried length r^. We chose eight values of ro, 
giving us D/r between 5 and 40 for the AEOS geometry, 
and created 100 independent realizations of Kolmogorov- 
spectrum phase screens for each. We then calculated cj> w 
over the AEOS aperture for each phase screen realization. 

Figure shows the expected trend of increasing vari- 
ance of (j) w with increasing D/r^. For cases where -D/ r o is 
greater than the number of WFS subapertures, the mean 
measured 4> w can be as high as (j) w = ±0.1 radians. More 
typical values are between ±0.04 radians. As we will 
shown in section 15. 1.21 waffle mode error starts to cause 
a significant reduction in Strehl ratio when phase vari- 
ations due to waffle grow beyond 0.1 radian. This sug- 
gests that WFRs which do not account for changing sky 
conditions throughout the night or which do not penal- 
ize waffle modes may not be able to accommodate these 
variations, and may allow waffle and waffle-like modes 
to build up, effectively enhancing the atmospheric waffle 
present in the observations. 

3.1.2. Effects of Waffle Mode on Strehl Ratio 

To analyze the intrinsic effects of waffle mode phase 
error on the observed PSF, we assumed perfect AO cor- 
rection (i.e., 4>ao — 0) with an introduced waffle er- 
ror of <f> w = 0.01 radians of phase difference. This 
corresponds to a typical value of <j) w present in our 
simulated Kolmogorov-spectrum phase screens with the 
AEOS WFS subaperture geometry. The result is a nearly 
perfect PSF attenuated by the regular pattern of satel- 
lite PSFs characteristic of fundamental waffle mode error. 
Here we note other waffle and wafflc-likc modes will pro- 
duce additional patterns of satellite PSFs in the image 
plane (M. Britton 2004, private communication). We 
note in passing that the strengths of these additional 
modes will depend on the details of the WFRs used. 

An example of our "waffle only" simulations can be 
seen in Figure where the waffle mode phase error has 
been increased to 0.5 radians to show the grid of satel- 
lite waffle PSFs. We note the primary, secondary, and 
tertiary waffle PSFs appearing at the intersection of grid 
lines a distance n\/2d (where n = 1, 3, 5, . . .). from the 
optical axis. As discussed in Section l2~3l the waffle PSFs 
which fall the grid points where n = 2, 4, 6, . . . hit the ze- 
ros of the sampled function sinc(dd), and do not appear 
in the final image. 

We define a waffle amplitude factor W to describe the 
enhancement of the waffle mode beyond the atmospheric 
4> w calculated from our simulated phase screens. In re- 
ality, W describes the amplification of the waffle mode 
due to the choice of wave front reconstructor, imperfect 
WFS flat fields, the effects of immobile DM actuators, 
and so on. 

We produced a series of PSFs with total phase errors 
ranging from 0.01 radian and 1.0 radian (peak-to- valley) 



by varying W between 1 and 100, keeping 4>ao = and 
<p w = 0.01. We then calculated the Strehl ratios of these 
PSFs by measuring the peak pixel value of the centroided 
PSF in the reconstructed aberrated wave fronts as com- 
pared to that from an image formed by a perfect wave 
front with the same image plane sampling. 

We present our Strehl ratio calculations in Table ^ 
We see little degradation in Strehl ratio for the cases 
with total waffle mode phase errors less than ~ 0.10 ra- 
dian - Strehl ratios in each of these examples remained 
higher than 0.99. Even for cases with 0.30 radians of 
total waffle, the calculated Strehl ratio remained above 
0.90. It is only when the waffle mode phase difference 
grows beyond this value that pure waffle mode error be- 
comes a major factor in reducing the Strehl ratio. This 
suggests that in the absence of other phase errors, fun- 
damental waffle mode error may not be an important 
factor in limiting the correction (and hence the observed 
Strehl ratio) in AO imaging. However, as more high- 
order AO (and ExAO) systems become prevalent, the 
use of WFRs that penalize or limit the growth of waf- 
fle modes in the AO system may be key to maintain- 
ing stable PSFs during AO observations. Additionally, 
the development of alternative methods of wave-front 
sens ing, such as the spatially filtered wave-front sen- 
sor llPov ne^r__^^dacinteB|2004|) or th e pyramid sensor 
(Raaazzoni 1996t lVerinaud et alJl2005]l . may be required 
to achieve the contrast levels necessary to achieve ExAO 
science. 

3.1.3. Symmetric and Antisymmetric Waffle PSF 

Terms 

In 12.3.11 w e presented the first two terms of the PSF 
expansion of iSivaramakrishnan et aLl l)2002|K and noted 
that the symmetric terms of this expansion would dom- 
inate the PSF as waffle mode phase errors increased. To 
examine this, we made use of the simulated PSFs de- 
scribed in 13.1.21 and differenced these PSFs with a per- 
fect PSF (po), generated with no phase errors. The result 
for three of these cases (4> w = 0.01, 0.03, and 0.10) can 
be seen in Figure [7| for aberrated PSFs with calculated 
Strehl ratios of 99.99%, 99.91%, and 99.0% respectively. 

We can see that power in the satellite waffle PSFs is 
accounted for by the P2, strehl term, 

p2,Strehl(®) = -\[a{a* *$**$*) + a* (a * $ * $)] (20) 

as illustrated by the on-axis "hole". There is some evi- 
dence of the antisymmetric p\ term in the diagonal pat- 
terns in these difference images, but their effect on the 
waffle PSF is small at these Strehl ratios. This anti- 
symmetry is more obvious as Strehl ratios decrease to 
between 97.0% and 90.0%, but the p\ term never domi- 
nates the intensity distribution. 

3.2. Simulating AEOS Adaptive Optics 

The AEOS WFS is a Shack-Hartmann system with 
an array of 84 /im lenslets, with each lenslet producing a 
Hartmann spot that is imaged onto a 4 x 4 group of pixels 
on the WFS CCD. The deformable mirror is a 31.5 cm 
diameter Xinetics DM supported by 941 lead magnesium 
niobate (PNM) actuators spaced 9 mm apart on a square 
grid, with 35 actuators spa nning the diameter of the DM 
^Roberts fc Nevmanl 12002]) . 
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Because of initial concerns about pupil wander during 
AO imaging, early AEOS wave front reconstructors did 
not use the entire pupil for wave front reconstruction. 
Subapertures near the edge of the pupil were assumed 
not to be fully illuminated during each observation, and 
were considered unreliable for active wave front recon- 
struction. To account for this, we modeled 33 active ac- 
tuators {N act = 33) across the pupil diameter instead of 
the full 35 actuators across the DM, giving us an under- 
sized DM with a 28.8 cm clear area and providing us with 
a projected actuator spacing of order 0.11 m/actuator at 
the primary mirror. 

To mimic finite DM actuator throw, we assumed an 
actuator stroke limit of ± 2.00 /im. While our sim- 
ulations model 855 active actuators on the primary 
(33 2 actuators x 7r/4), the actual AEOS WFRs actively 
control 810 actuators, with actuators projected near the 
edges of the primary and secondary mirrors slaved to 
neighboring actuators. We did not model slaved actua- 
tors, nor did we model the effects of actuator influence 
functions, non-linear actuator motions, or actuator hys- 
teresis. We also do not account for system response time 
or calibration errors. 

3.2.1. Seeing Measurements at Haleakala 

The AEOS Adaptive Optics imaging data we used for 
our comparison comprises four sets of twenty-five closed- 
loop /-band observations of the K3 II star HR 7525 (7 
Aql, V — 2.72) taken at two separate exposure times (48 
ms and 73 ms) with the Visible Imager (Vislm) camera 
using a 2 magnitude neutral density filter. This star 
also served as the WFS target during these observations. 
These data were acquired on 8 October 2002 and span a 
total of 5.43 minutes from the first to the last exposure 
in the sequence. The images were processed (dark, bias, 
and flat field corrected) using the task CCDPROC in 
IRAF. 5 

Measurements of ro taken with the University of 
Hawaii's Day Night Seeing Monitor (DNSM) during the 
night of our observations showed a mean measured r$ = 
11.44 ± 3.54 cm at A = 0.500^m from 320 data points 
over 9.22 hours. The median measured ro for this night 
was 11.97 cm. The mean zenith-corrected ro for this 
night was ro = 13.17 ± 3.88 cm, with a median value of 
13.73 cm. The nine measurements acquired just before, 
during, and after our observations showed a mean zenith- 
corrected ro of 11.63 cm with minimum and a maximum 
measured values of 8.84 cm and 13.60 cm respectively. 

This variation is apparent in our data, as measure- 
ments of the full- width at half maximum (FWHM) of the 
PSF in each image showed a median of 3.17 ± 0.26 pixels 
for 92 images, with a minimum FWHM of 2.77 pixels 
and a maximum FWHM of 4.49 pixels. Over its nominal 
optical bandpass, (from 0.70 fim to 1.050 fim) the Vislm 
has a pixel scale of 0.0211 ± 0.004 arcsec/pixel, giving us 
a median AO FWHM of 0.0669 arcsec. 

Unless otherwise noted, we have adopted the value 
ro = 11.60 cm to define our nominal seeing conditions. 
We scale the value of ro J-band and to //-band following 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



the relationship 

/A \ 6/5 

r i= r o\Y ) cos3/5c ( 21 ) 

where we assume the zenith angle, C = l|Hardvl fl998). 
In this way, we calculate r\ = 22.92 cm and rn = 
48.20 cm, and use these values in our simulations. 

3.2.2. Estimating Waffle in AEOS I-band Data 

To estimate the amount of waffle present in the AEOS 
/-band data, we first required a knowledge of the Strehl 
ratios of the images themselves. We measured Strehl 
ratios by first simulating the diffraction pattern of the 
AEOS telescope, as computed from the analytical ex- 
pression for the diffraction pattern of a centrally obscured 
circular aperture. We extracted photometry of the PSF 
and diffraction pattern, and registered the PSF and the 
diffraction pattern to sub-pixel resolution by cross cor- 
relation. The PSF was then shifted to have the same 
sub-pixel center as the diffraction pattern, and the Strehl 
ratio is computed f rom the ratio of peak pixel intensities 
in both images (see lRoberts et all ()2004|) for a discussion 
of Strehl computation algorithms). Strehl ratios mea- 
sured in this manner ranged from 0.17 to 0.28. 

We measured the intensity of the peak pixel in each of 
the four satellite "waffle PSFs" observed in the images 
of HR 7525, and compared these measurements to the 
peak pixel of the on-axis PSF. We found that the peak 
intensity of each waffle PSF was of order 0.25% the peak 
intensity of the on-axis PSF. The total flux measured in 
the waffle PSFs was ~ 1.6% the flux of the on-axis PSF. 

We then compared this result with the relative inten- 
sities of the waffle spots produced through simulation. 
Our simulations made use of five independent realiza- 
tions of Kolmogorov-spectrum phase screens, each with 
ro = 22.92 cm. Our AO correction model used the 0.9 
power law filter in spatial frequency space, as described 
earlier in Section |31 The residual phase that remained 
after AO correction {i.e., <Pao 7^ 0) was then used to sim- 
ulate a PSF for each of the five incoming phase screens. 
These PSFs where then co-added to produce our final 
image. 

We simulated a set of monochromatic PSFs using this 
set of five phase screens, calculating <f> w from the phase 
screens themselves. We chose waffle amplitudes W rang- 
ing from zero, to atmospheric, to the extreme {e.g., 
100W), and ran simulations using the same five input 
phase screens for each value of W. We then examined the 
resultant images, and assumed the peak intensity of the 
simulated waffle PSFs corresponded to the peak wave- 
length of the /-band filter transmission function used in 
the Vislm observations. We then chose the value of W 
which produced differences in the peak waffle spot and 
central PSF intensities which best matched the observed 
values from the data. 

Our simulated images are generated at a pixel scale of 
A/4D, or 0.0125 arcsec/pixel, nearly a factor of two finer 
sampled than the true Vislm sampling at this wavelength 
(0.0211 ± 0.004 arcsec/pixel). We rebinned the simu- 
lated data to matc h the Vislm pixe l scaling using the 
IDL task FREBIN <|Landsmanl 11993). We added Gaus- 
sian detector read noise (12 rms) and dark current 
(22 e~ pixel -1 s _1 ) into our rebinned PSFs, and mod- 
eled Vislm detector charge diffusion and optical effects 
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by convolving our simulated images with a Gaussian pro- 
file of a = 1.15 pixels using the IRAF task GAUSS. This 
has the effect of limiting the maximum observable Strehl 
ratio in simulated Vislm data to S ~ 0.70, which matches 
expectations for the Vislm and its optics as used in Oc- 
tober 2002. 

We found a waffle mode amplitude of 16W at A = 
0.88 /im best matched the intensity of waffle spots 
in these Vislm /-band images. An example of this 
monochromatic simulation with full AO correction, both 
without and with waffle, can be seen in Figure |S1 We 
used this as the starting point to extend our monochro- 
matic simulations to full 7-band simulations. 

3.2.3. Simulating I-band Images 

We simulate a set of J-band images by first generat- 
ing thirty-six monochromatic images across the band- 
pass, from A = 0.700 /im to 1.050 /Ltm (the effective cut- 
off wavelength of the Vislm detector) in increments of 
0.01 /xm. We use the same Kolmogorov-spectrum phase 
screens for each image at a given A;, scaling the phase 
difference of each phase screen relative to the bandpass 
central wavelength e.g., 4> a tm • (Aj/A re /), where A re / is 
the peak wavelength of the filter transmission function 
in the bandpass. We maintained the same 0.9 power 
law AO filter in spatial frequency space as before, with 
N ac t — 33 actuators across the primary and a waffle 
mode amplitude of 16W. We did not consider the effects 
of slaved or dead actuators as part of our simulations, 
but did include a ±2.0 /im DM stroke limitation in our 
analysis. 

We rebinned the pixel scale of each monochromatic 
PSF by the ratio of the simulated wavelength to the 
reference wavelength to account for the wavelength- 
dependent resolution changes inherent in Fourier meth- 
ods. We then multiplied each rebinned \ PSF by the 
transmission expected for a standard Johnson/Bessell 
/-band filter transmission function obtained from SYN- 
PHOT in STSDAS, 6 and combined the set of monochro- 
matic PSFs for a given phase screen into a single broad- 
band PSF using the IRAF task IMCOMBINE. 

We generated five broadband PSFs in this manner, and 
co-added these PSFs into a single "short exposure PSF." 
Assuming a speckle correlation time of order a few tens 
of milliseconds, these five PSFs, when averaged, would 
approximate a ~ 100msec PSF. We note that longer 
exposures can be simulated in this way by adding sta- 
tistically independent realizations of PSFs together, but 
would lack any correlated or longer-term effects (such 
as quasi-static speckles) without the addition of static 
sources of phase errors. 

With the methods described above, our our simulated 
Kolmogorov wavefronts produce PSFs with Strehl ra- 
tios between 0.544 to 0.635 with purely atmospheric 
phase disturbances. When waffle mode phase error with 
W = 16 is introduced, the Strehl ratios decrease to be- 
tween 0.406 and 0.556 for the same phase screens. When 
convolved with the Gaussian filter, the Strehl ratios of 
our simulated PSFs were further reduced by 70%, to 
between 0.284 and 0.389. However, the resultant final 

6 SYNPHOT and STSDAS are products of the Space Telescope 
Science Institute, which is operated by AURA under contract to 
the National Aeronautics and Space Administration. 



images produced waffle spots with the same morphology 
and intensity distribution as in Vislm images (Figure EJ- 
The total flux in the visible waffle spots, ~ 1.8% the flux 
of the central PSF, is in good agreement with the ~ 1.6% 
value measured in the Vislm data. The peak intensities 
of each waffle spot, ~ 0.25% of the peak of the central 
PSF, is also in agreement with the Vislm data. 

While the simulated Strehl ratios remain higher than 
expected for the ro values observed in the Vislm im- 
ages from this night, we note that our simulations only 
account for two sources of phase error: residual atmo- 
spheric variations and the fundamental waffle mode. We 
do not account for pointing variations or pupil wander in 
the telescop e beam, nor do we employ a realistic model 
of the DM. lAbreu et all (|2000j) noted that pinned and 
slaved actuators may provide as much as a 10% hit to the 
observed Strehl, while actuator influence functions and 
actuator dynamics will also limit the quality of the wave- 
front reconstruction. Any observed non-common path 
aberrations between the AEOS WFS/DM and the Vislm 
will also degrade the observed Strehl ratios, as will the 
effects of the Vislm detector pixel modulation transfer 
function (MTF). 

We stress here that these Vislm images, acquired early 
in the AEOS AO system lifetime, did not benefit from the 
waffle-suppressing WFRs that are now in use at AEOS 
(G. Smith 2003, private communication). Images ob- 
tained more recently in I— , J—, H-, and K— bands have 
shown waffle mode error in AEOS imaging to be at or 
near nominal (atmospheric) levels l|Perrin et al.ll2004t) . 

3.2.4. The Magnitude of Waffle Mode Phase Error 

Following iSandler et alJ l)1994l) . we note that for 
Strehls > 0.30 the final observed Strehl ratio of a sys- 
tem can be considered as a multiplicative combination 
of Strehl ratios due to successive uncorrelated sources of 
error. Thus 

S = Sfr... = e-<i e-< . . . e -»3„ (22) 

where is due to one class or family of aberrations 
(i.e., high-order aberrations), while a^ 2 is due to another 
class of aberrations uncorrelated with the first, etc. 

Neglecting the optical effects of the Vislm, we can cal- 
culate the mean square wavefront error for atmospheric 
phase disturbances, <r? Q , from the mean Strehl ratio of 
our simulated images, 

S = S,p AO = er a *Ao = 0.589. (23) 

If we then assume a^, AO = 2irz/\, where A = 0.880/im, 
we find z - the residual RMS wavefront of the AO- 
corrected atmospheric phase screen - to be 101.9 nm. 
If we then consider the addition of waffle mode phase 
error with W = 16, we find, 

S = S+ AO S4 w = e^oe-^. = 0.488. (24) 

With S<j, AO = 0.589 as above, = 0.829, and thus 
cr^ = 0.188. Following the method above, we calculate 
the residual RMS wavefront due to waffle mode phase 
error to be z w = 60.7 nm. 
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3.3. Other Potential Source of Waffle Mode Phase 
Errors 

While our simulations were able to reproduce the gen- 
eral morphology of the early AEOS AO system PSF as 
well as the shape and energy seen in the primary waffle 
PSFs, we realize our work does not fully reproduce the 
AEOS optical system, WFS geometry or WFRs. In fact, 
the presence of immobile DM actuators in the AEOS 
pupil imprin ts a phase and an a mplitude error across the 
entire pupil. lAbreu et alJ 1)20001) has estimated the Strehl 
hit to AEOS AO system performance due to these im- 
mobile actuators at ~ 10%. Initial work has suggested 
these failed actuators can mimic localized waffle modes 
described bv lC^ell l(200l . 

Given a DM surface map taken with an interferome- 
ter, it is relatively straightforward to calculate the waf- 
fle present at the DM surface. Comparing this to the 
strength of the waffle seen in real images will indicate 
whether WFRs tailored to minimizing the effects of failed 
actuators will help to significantly improve the image 
quality. A reconstructor which drives the wave front 
to the best-fitting plane through the failed actuators' 
heights on the DM is one obvious strategy. However, 
such a strategy would adversely affect the usable stroke 
of the DM, and thus would only be beneficial only when 
the seeing is good. 

Circular aperture boundaries and partial subaperture 
obscuration by edges and secondary support spiders will 
affect the equations relating WFS slopes to actuator com- 
mands. Other waffle modes can be introduced by sec- 
ondary support structure obscurations. In the extreme 
case, a completely obscured line of subapertures parti- 
tioning the pupil into disconnected sections can result 
in free-floating piston between the sections in the re- 
constructed wave front. While this form of error might 
be less obvious, it can adversely affect Strehl ratios and 
dynamic range of companion searches and debris disk 
studies because of its proximity to the core of the re- 
constructed PSF. These quadrants must somehow be 
stitched together explicitly by the reconstructor alge- 
bra with the mathematical equivalent of extra equations 
of constraint. Partially- illuminated subaperture data 
should be ignored or used differently from WFS centroid 
data from fully-illuminated apertures in order to reduce 
their contribution to undesirable behavior of the recon- 
structor. 

4. PREDICTING AEOS AO SYSTEM PERFORMANCE FOR 
NEAR-INFRARED OBSERVING 

In order to understand the potential capabilities of 
diffraction-limited near-infrared coronagraphy on AEOS, 
we extrapolated our /-band simulations to //-band em- 
ploying the same power-law filter in spatial frequency 
space to simulate wave front correction as we did for 
our /-band simulations. The use of the same AO repre- 
sentation provides a reasonably accurate simulation of a 
facility AO system working independently of the science 
imaging instrument, as in the case of the AEOS AO sys- 
tem (among many others). Here, we discuss details of 
our //-band simulations. 

We generated realizations of Kolmogorov-spectrum 
phase screens with ro = 48.20 cm, which created a set of 
//-band phase screens with an effective D/ro of 7.5. We 
assumed a waffle amplitude of 16W as before, with a DM 



stroke limit of ± 2.00 /im, and ran a set of monochromatic 
simulations covering the wavelength range from 1.43 fjxa 
to 1.83 /im in increments of 0.01 /im. We combined this 
set of monochromatic images into a single broadband im- 
age (with the appropriate rebinning, described in l3.2.3|l . 
weighting each wavelength's PSF with the Mauna Kea 
Obse rvatories (MKO) //-band filter transmission func- 
tion ijSimons fc Tokunagal2002D . We did not include any 
other optical characteristics of the telescope, AO system, 
or imaging system for these simulations. 

Our initial simulations suggested the AEOS AO system 
is capable of high Strehl ra tio images in the //-b and even 
without improved WFRs (|Makidon et al]l2003ft . though 
the presence of significant satellite waffle spots could 
make faint companion searches problematic in particu- 
lar regions relative to a bright AO target star (see Fig- 
ure However, recent improvements to the AEOS 
WFRs (G. Smith 2003, private communication) have re- 
duced the waffle mode error remaining in AEOS imaging 
to levels that are almost undetectable in typical images. 

In March 2003, J— , H— , and /C-band images of several 
stars were taken with the Kermit Infrared Camera as 
part of a study of atmospheric turbulence. This data 
presented us an opportunity to validate our simulations 
against act ual data. The cam era and data set up are 
detailed in IPerrin et al.l l|2004|) . Our simulations with 
a waffle amplitude of 16W proved to be too aggressive 
for the AEOS WFRs in use at the time of the Kermit 
observations, and we have since re-determined our "best" 
W for AEOS //-band imaging to be closer to 4W, well 
into the noise of most imaging systems. 

On-the-sky results using the Kermit Infrared Cam- 
era have since demonstrated the improvements to the 
AO system, showing only the faintest traces of waffle 
mode phase error in over-exposed images. In Figure ITT1 
we show Kermit short- //-band data of the G8III star 
HR 4667 (V = 4.94) acquired at AEOS on 16 April 2003. 
Images shown were taken with two exposure times: 1.5 
seconds (at left), and 5.0 seconds (middle), with mea- 
sured Strehl ratios of order ~ 0.50. Here, we found that 
Kolmogorov-spectrum phase screens with ro = 24.2 cm 
at //-band best matched the general morphology and 
Strehl ratios of these images, suggesting the seeing on 
the night these data were acquired was below nominal. 
However, we note that these relatively long exposures 
will undoubtedly contain many tens of speckle lifetimes, 
while our simulated images only use five independent re- 
alizations of a simple atmosphere. As such, atmospheric 
speckles will tend to be "washed out" in the real data, 
b ut instrumental effe cts will still be visible. 

IPerrin et al.l l|2004|) reported measured //-band Strehl 
ratios as high as 83% during good seeing conditions at 
AEOS. Our simulations match these Strehl ratios in the 
absence of other aberrations, though are yet unable to 
reproduce the quasi-static speckle structure present in 
the AEOS /-band and //-band images. In the future we 
hope to be able to determine phase and amplitude maps 
for the AEOS AO system beam as it is presented to the 
Lyot Project coronagraph and to the Kermit Infrared 
Camera to enable proper PSF calibration both with and 
without the coronagraphic occulting and Lyot stops in 
place. 
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5. CONCLUSION 

One of the goals of this work is to show the util- 
ity of simulating AO system performance in a simple 
manner, rather than employing detailed, component-by- 
component modeling, to enable long exposures to be sim- 
ulated in a short time using modest computing power. 
We presented the AEOS telescope and its AO system as 
a specific case of interest, but note this method should 
work on almost any system given imaging data with 
which to calibrate the spatial filter model of AO. Other 
software tools, such as PAOLA ( Performance of A dap- 
tive Optics for Large Apertures, Jolissant (2004)), in- 
corporate more physics into such "fast" methods of AO 
system simulation. These types of simulations are useful 
both for highlighting areas where the AO system perfor- 
mance could be improved, and for predicting instrument 
performance to investigate observing scenarios to inves- 
tigate what science is within the scope of an AO system 
and its science camera. 

AEOS and its AO system can prove to be an excellent 
testbed for technologies requiring diffraction-limited ob- 
serving. A better understanding of the AO system, such 
as that attempted by this analysis, will allow the as- 
tronomical community to build instruments to take bet- 
ter advantage of the unique capabilities of AEOS. Stud- 
ies like these will provide the astronomical community 
an opportunity to learn lessons which can be applied to 
high-order AO systems on 8 and 10 m class telescopes. 
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Fig. 1. — The unit cell of the Fried geometry wave front sensor. Actuators are shown as solid squares, and lcnslct edges as dashed lines. 
The (x, y) coordinate system is in wave front sensor pixels. The centroid of the image produced by the lenslet is shown as the dark circle. 
The subaperture, when projected back on to the primary, is a square of side a. 




Fig. 2. — A square aperture Fried geometry wave front sensor array, with sixteen actuators (shown as squares, each at a height hi) 
and nine wave front sensors. Wave front slopes Si = (s x i,s y i) denote the instantaneous centroids of a lenslet spot in the Shack-Hartmann 
subaperture (in on-chip pixels). 
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Fig. 3. — The AEOS pupil and subapcrturcs over which we calculate our simulations. Left: the telescope input pupil. The diameter of 
the telescope is 512 pixels, with the secondary mirror obstruction subtending 56 pixels. Right: the distribution of subapertures over the 
telescope pupil over which we calculate the mean 4> w in the incoming phase screen. 




Fig. 4. — Simulated Kolmogorov-spectrum and AO-corrected phase screens over the telescope pupil. Left: the incoming phase screen 
over a portion of the AEOS pupil with linear scaling. The AEOS secondary obstruction is shown at the upper-right. Center: the residual 
AO-corrected phase screen over the same portion of the AEOS pupil. Here, the linear scaling has been enhanced relative to the figure at 
left to show the residual structure. Right: the same residual AO-corrected phase screen as at center with waffle mode introduced, shown 
at the same linear scaling as at center. Here, <j> w = 0.03 radians with W = 10. 
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Fig. 5. — Measured mean atmospheric if>w values from 100 Kolmogorov-spectrum phase screens with the AEOS telescope and wavefront 
sensor subaperture geometries (top), with the mean measured value and standard deviation of the mean (bottom). 





Fig. 6. — Simulated monochromatic image assuming the AEOS telescope and WFS geometries with fundamental waffle mode error 
defined by <j> w = 0.01 radians and W = 50. Atmospheric correction is assumed to be perfect {i.e., no residual atmospheric phase error 
other than waffle). This image is on a min-max scaling with a logarithmic stretch. In this extreme case, one can clearly see the primary, 
secondary, and tertiary waffle spots around the central PSF (see text). 
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Fig. 7. — Difference between simulated AEOS monochromatic PSFs with and without waffle mode error. From left to right: PSFs with 
waffle 4> w = O.Of, 0.03, and 0.10 can be seen in the PSFs from left to right with calculated Strehl ratios of 99.99%, 99.91%, and 99.0% 
respectively. In each case, atmospheric correction is assumed to be perfect (i.e., no residual atmospheric phase error other than waffle). 
Subtracting the perfect PSF results in a "hole" where the on-axis PSF with waffle would be, indicative of a reduced Strehl ratio in the 
waffle images. The antisymmetric component of the PSF due to waffle mode error is visible, though the symmetric component of the waffle 
mode error dominates. 




Fig. 8. — Simulated AEOS monochromatic PSFs without and with waffle mode error, displayed on the same logarithmic scaling. These 
images assume D/tq = 32 at A = 0.88 /im with a power law AO correction with exponent 0.9 (see text). Neither read noise nor Poisson 
noise arc included here. 



* 




♦ 





















Fig. 9. — At left, an AEOS Adaptive Optics /-band image of HR7525 acquired with the Visible Imager (Vislm) on 08 October 2002 (UT). 
Strehl ratio of this image has been estimated to be ~ 20%. At right, a simulated AEOS Adaptive Optics /-band image assuming waffle 
mode error. The Strehl ratio of the simulated image was calculated to be ~ 60% before convolution with a Gaussian filter and and the 
addition of Poisson noise, which reduced the calculated Strehl to ~ 42%. Other features in the observed PSF can be explained by the 
addition of phase errors in the input pupil. 
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Fig. 10. — Simulated AEOS ff-band PSFs. Both of the PSFs shown here were co-added from five individual PSFs, each using a separate 
realization of Kolmogorov-spectrum phase screens with ro = 45.6 cm as input. At left is a PSF with W = 16, with W = 4 for the PSF 
at right. Poisson noise and detector read noise have been added to these PSFs. Calculated Strehl ratios are 0.79 for the PSF at left, and 
0.82 for the PSF at right. 




Fig. 11.— Kermit Infrared Camera short-Jf-band images of the G8III star HR 4667 (HD 106714, HIP 59847) acquired on 16 April 2003 
with AEOS Adaptive Optics. At left, a single 1.5 second exposure with a measured Strehl of ~ 0.51. A 5.0 second exposure of the same 
star is shown at middle. At right is a simulated PSF with a single-layer Kolmogorov-spectrum phase screen with D/ro = 15.0, with a 
calculated Strehl ratio of ~ 0.55. Since our simulations assume perfect optical surfaces, the calculated Strehl ratios for the simulated PSF 
are higher than those imaged with Kermit on the night of the observations. 



TABLE 1 

Average Strehl Ratio for <j> w = 0.01 Radians 



Waffle (W) 


<Ptotal (rad) 


Strchl a (%) 


Waffle (W) 


<t>total (rad) 


Strchl a (%) 


0.0 


0.00 


1.0000 


20.0 


0.20 


0.9605 


1.0 


0.01 


0.9999 


25.0 


0.25 


0.9388 


3.0 


0.03 


0.9991 


30.0 


0.30 


0.9127 


5.0 


0.05 


0.9975 


50.0 


0.50 


0.7702 


10.0 


0.10 


0.9900 


75.0 


0.75 


0.5354 


15.0 


0.15 


0.9776 


100.0 


1.00 


0.2919 



a Strehl ratios in TableHlwere computed by taking the ratio of the peak intensity in the 
aberrated image with that of a perfect image. 



